Stimulated Raman Transfer of SIT Solitons 



B.D. Clader and J.H. Eberly 
Department of Physics and Astronomy 
University of Rochester 
Rochester, NY 14627 



Abstract 

We treat theoretically the two-pulse plane-wave propagation prob- 
lem associated with the stimulated Raman process. A strong pump 
pulse and a weak probe pulse interact self-consistently during mutual 
propagation through a medium in the three-level lambda configura- 
tion. We treat the pulse interaction fully coherently, but take initially 
incoherent mixed ground states to describe the medium. We present 
analytic pulse envelope solutions to the nonlinear evolution equations, 
and then test them numerically. We show that the entire process can 
be interpreted consistently as a three-regime evolution, in which the 
first and last regimes display behavior quantitatively the same as self- 
induced transparency. They are connected through a Raman transfer 
process in the intermediate regime. 



1 Introduction 

It is a pleasure to join the other authors of this special volume in celebrating 
the exceptionally wide variety of accomplishments in the distinguished career 
of Professor Girish S. Agarwal. His impressive contributions as a theoretical 
physicist in quantum optics began with more than two dozen papers written 
in our Department in Rochester while still a student pj. It seems certain 
we will continue to benefit from many further interesting calculations and 
valuable insights in years to come. 

As coherence in quantum optical contexts has provided a constant source 
of fruitful inspiration to Professor Agarwal, we've chosen the Raman effect 
during fully coherent pulse propagation in three-level media as the topic of 
our contribution. This is an area that has prompted recent studies of his 
own [2 El m El m [7] , and is an area of wide theoretical and experimental 
activity at present. The door to this domain was opened by the famous work 
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of McCall and Hahn [8j on two-level effects in ruby, in which they found 
theoretically and demonstrated experimentally the first optical soliton, the 
27r sech pulse of self- induced transparency (SIT). 

The coherent pulse propagation domain was gradually expanded to in- 
clude pulse-pair propagation including three- level solitonic behavior pilUlfTUl 
[TT\ [12], and the phenomenon of electromagnetically induced transparency 
(EIT) then emerged in the seminal work of Harris [13]. This has led to 
exciting discoveries including superluminal and ultra-slow pulse propaga- 
tion. A consistent theoretical challenge, particularly in numerical modeling 
of short pulses, has been the effect of inhomogeneous broadening (see some 
treatments in [BJ [Qj [HI [E]), and that is one element of the present paper. 

The Raman effect was the first effect in physics of a specifically three- 
level character, and it is the basic model for essentially all experiments on 
lambda systems with two optical fields. Our contribution can be considered 
as an examination of the cooperation between SIT and the Raman process, in 
contrast to the study by Agarwal, Dey and Gauthier ^ on the competition 
between EIT and Raman processes. Of course the stimulated Raman effect 
(STR) has been actively studied for decades (for an early review see |16j). 
The STR is a process by which a strong laser beam incident on the "pump" 
transition of a three-level lambda system causes amplification of a weak 
"probe" pulse injected at the Stokes frequency (see Fig. 1). Traditionally 
STR has been studied in atomic and molecular systems, but it has also been 
examined in Bose condensates [3j. 

In this paper we present analytic and numerical solutions to the coupled 
plane- wave Maxwell-Bloch system for stimulated Raman scattering, without 
the usual assumption of undepleted pump pulse and with careful considera- 
tion of non-trivial changes in both the probe and pump pulse shapes due to 
propagation. We show in the fully coherent case (pulses short enough to ne- 
glect homogeneous damping) that complicated exact analytical results can 
be interpreted relatively easily, partly in SIT terms, by introducing three 
"regimes" of interaction. We show how the three regimes can be considered 
separately despite the tightly coupled nonlinear nature of the solutions. 

2 Theoretical Raman Model 

We consider the coherent propagation of two laser pulses in the x direction 
in a medium of three level atoms in the lambda configuration as shown in 
Fig. [H The electric field of the individual laser pulses can be written as 
Ea{x, t) = £a{x, t)e~*('^'»^~'^'»*) + c.c, where £a is the slowly varying envelope 
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of the electric field, ka is the wavenumber, and uja is the frequency (similarly 
fovEb{x,t)). 
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Figure 1: Three-level atom in two-photon resonance with level 1 connected to level 
3 via the Rabi frequency fla and level 2 connected to level 3 via the Rabi frequency 
fib- Each laser field is detuned from the excited state resonance by an equal amount 
A, and the excited state decay rate is given by 73. 

The Hamiltonian of the system in the rotating wave picture is given by 



where h.c. refers to the hermitian conjugate of the previous operator. For 
the 1-3 transition we have = 2da£a/^, da is the dipole moment, and 
A = — (wi + cja), and corresponding notation for the 2-3 transition. 
Note that the same A is the one-photon detuning of each transition below 
resonance. 

We will describe each individual atom by its density matrix. The von 
Neumann equation for the density matrix ifi^ = [H, p] gives the following 
density matrix equations, with the use of the Hamiltonian in Eq. ([T|): 





pn = i^Psi - i^Pn 



(2a) 



.^b .^b 
P22 = «^P32 - i-:^P23 




(2b) 



(2c) 



P12 = i^P32 - i^PlS 



(2d) 




P23 = ^-^{P33 - P22) - 'i-^P21 + ^Ap23- 



(2e) 



(2f) 
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Eqns. ([2|) assume that the duration of the laser pulses is sufficiently short 
that the excited state decay rate 73 can be neglected because we consider 
pulses with durations much shorter than I/73. Maxwell's equation for the 
field gives these two slowly varying envelope equations 



where /x^ = NdatOa/^^o is the atom field coupling parameter, is the 
density of atoms, da is the dipole moment matrix element between atomic 
states 1 — > 3, and similarly fih for states 2 — > 3. The notation (...) = 
J^(. . .)F{A)dA accounts for inhomogeneous broadening by averaging over 
a range of detunings A with the function F(A) = (r2*/\/27r)e~*^"^2 ^)V2^ 
where T| is the inhomogenous lifetime. For these equations to permit an- 
alytic solutions, one needs /i^ = /i;, = /i, which we will assume for the 
remainder of this paper. 

3 Pulse Propagation Analytic Solutions 

We solve Eqns. and ([3]) using a method developed by Park and Shin 
|12j which is based on the Backlund transformation. We assume the initial 
density matrix of each atom is in an incoherent mixture of the two ground 
states, which is written in explicit form as 




(3b) 



(3a) 




a 













The solution to the field envelopes is given by 
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(5a) 




(5b) 
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where r is the nominal pulse width and T = t — x/c \s the retarded time, 
and A and B are given by 

i? = l/3|^^^, (6b) 

where the function g{T, ) is the Beers absorption coefficient (inverse Beers 
length) sometimes written a/) in the case of Doppler broadening: 

^ r FjAJdA _ 



Related solutions were previously given in [10] and [12] for the case where 
|ap = 1, and in the limit where inhomogeneous broadening is ignored. 
For the remainder of this paper we will be considering the opposite limit, 
corresponding to large inhomogeneous broadening as in the familiar SIT 
scenario [8], and a~^^ will be the appropriate length unit for spatial pulse 
evolution. 



4 Three-Regime Behavior 

The analytic pulse solutions are plotted in Fig. [2j Their solution formulas 
dlj) resist easy interpretation. However, they can be unravelled as follows. 
From the plots we see three distinct types of evolution. When only is 
significantly present we will speak of regime I (frames 1 and 2 of Fig. [2|). 
During "transfer" from pump to Stokes when both pulses are substantial 
we have regime II (frames 3-5). During regime III only Qb is significantly 
present (frame 6). We will isolate these three regimes analytically by looking 
at the asymptotic behavior of the pulse solutions. By considering x <^ —c^i) 
we can look at regime I, and by considering x ^ we can look at regime 
III. The transfer regime II corresponds to x ~ 0. 
The Area of a pulse is defined to be 

/oo 
n{x,t)dt, (8) 
-oo 

from which the pulse Areas of solutions ([5|) can be shown to be 

0Jx) = --^rtan"^(/i(x)) +cot"^(/i(x))l (9a) 
h[x) 

0bix) = —^[tan-\hi-x)) + cot-i(/i(-x))] (9b) 
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Figure 2: Plots of the analytic pulse solutions given in Eq. The horizontal 

axis is X in units of the absorption length aj^^, and the vertical axis is the pulse 
Rabi frequency in units of t~^. The background is slightly shaded to indicate the 
presence of the lambda medium. The solid curve is ila and the dashed curve is 
rii,. The plot shows Raman amplification as an exchange process between two SIT 
solitons. The input pulse on the 1-3 transition amplifies a weak probe pulse on 



the 2-3 transition. Parameters: lal — 0.8, 



= 0.2, r « 3Tn* 



where 

/,(x) = (l + e(^-^)^y/2^ (10) 

In the asymptotic Hmit x <C —ct]^^ it is straightforward to show that 9a{x) ~ 
2ir and ^^(x) 0. In the opposite hmit, where x » cud^^ we have the 
opposite result, 9aix) ~ and Oi,{x) In. 

In regime I, Eq. ^ can be approximated by (recalling the definitions of 
A, B, and T): 

Q,a ^ — sech I — — \a\^—^x] (11a) 

T \ T 2 J 

~ 0, (lib) 

which corresponds to a single 27r sech-shaped pulse traveling at a group ve- 
locity ofvg = c/ (I+Act). In regime III the pulse solutions are approximated 
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by 



Qa^O (12a) 
sech - m'^x^ , (12b) 

which corresponds to a single 2tt sech-shaped pulse traveling at the group 
velocity of Vg = c/(l + Bct). Thus in both regimes I and III we see that 
the solution reduces to that of a 27r Area sech shaped pulse moving in a 
two-level medium. These asymptotic solutions correspond exactly to SIT 
solitons [8j. These two regimes are connected via the intermediate regime 
II where the three level character of the medium causes the Raman transfer 
process to occur, and the full solutions given in Eq. ([5]) must be used. 



5 Numerical Solutions 

To test the potential experimental relevance of the analytic solutions pre- 
sented, we have to replace the infinite uniform medium assumed to this 
point by a more realistic medium with finite edges. We then solve Eqns. ([2]) 
and dSI) numerically. The analytic solution shapes and areas should not be 
assumed in advance, so at the entrance face to the medium we will inject 
gaussian input pulses given by 

r^/2^ 

^t^ = -^e-^, (13b) 
rv27r 

where the parameters 6a and 9i, give the input pulse areas. In the previous 
section we showed how the analytic solution can be separately interpreted 
in three distinct regimes by considering asymptotic solutions. We now chose 
input pulses for the numerical solution specifically to test the validity of our 
hypothesis that this process can be broken up into the three regimes defined 
above. 

In Fig. [3] we plot the numerical solutions to Eqs. ([2]) and ([3]). We 
assume the medium is initially configured just as in Fig. [2l with |ap = 0.8 
and = 0.2, and that the medium is inhomogeneously broadening such 
that T 3T2*. We choose input pulse Areas 9a = 1.37r and 9b = O.OOSvr. 
The numerical solution shows that the 1.37r gaussian-shaped input pulse is 
reshaped into a 27r sech-shaped pulse in regime I just as a pure two level SIT 



7 



4 







5 ■ 















-40 -30 -20 -10 10 20 30 40 -40 -30 -20 -10 10 20 30 40 



Figure 3: Plots of the numerical pulse solutions of Eqs. ([2]) and ([3]). The horizontal 
axis is X in units of the absorption length ajj^ , and the vertical axis is the pulse Rabi 
frequency in units of t~^. The shaded zone indicates the location of the medium. 
The solid curve is and the dashed curve is fif,. The plot shows three distinct 
processes. A two level SIT process where the input I.Stt gaussian pulse is reshaped 
into a 2tt sech shaped pulse (frames 1-2), followed by a Raman amplification process 
where the input pulse on the 1-3 transition amplifies a weak probe pulse on the 2 - 
3 transition (frames 3-5), and finally a second two level SIT process where the weak 
probe pulse has become a 2tt sech shaped pulse (frame 6). Parameters: |ap = 0.8, 
= 0.2, r « 3T2*, 0a = 1.37r, Ob = O.OOStt. 



description would suggest. The Raman transfer process then takes place in 
regime II, which amplifies the weak O.OOSvr Stokes pulse to a 27r sech shaped 
pulse on the 2 to 3 transition. This final pulse in regime III is an SIT soliton, 
and will propagate without any dispersion or further modification. We also 
plot the numerically integrated pulse Area in the bottom frame of Fig. HI 
One can clearly see the input pulse as it is reshaped into a sech shape is also 
changed to 27r total Area, just as predicted by standard SIT theory. 

The Area of the input pulses used in the numerical solutions is an ex- 
perimentally controllable parameter. This value determines how many ab- 
sorption lengths into the medium the pulses will travel in regime I before 
beginnning the transfer process in regime II. The top frame of Fig. H] shows 
the pulse Areas given in Eqns. ([9]). The input area of the weak pulse de- 
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termines how far to the left of the asymptotes in Fig. H] we start. In the 
numerical solution we choose an input pulse area for pulse b of 0^ = O.OOSvr 
which corresponds to a^x ~ —20 in the analytic Area formula given in Eq. 
([9]). Thus within about 20 absorption lengths one would estimate the two 
pulses would be of roughly equal Area. In Fig. Owe indeed see the pulses in 
frame 3 are almost of equal Area, after traveling about 20 absorption lengths 
through the medium. This can be more clearly seen in the plots of the pulse 
Areas in Fig. [J] where we see the pulses have equal Area near aux ~ 20 
as predicted. Thus we can use Eqns. ([9]) along with the input Area of the 
Stokes pulse to determine where the Raman transfer will occur even when 
the input pulses are of different shapes. 
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Figure 4: Plots of the pulse Areas from Eqs. ([9]) in the top frame, and from inte- 
grated pulse Area from the numerical solution in the bottom frame. The horizontal 
axis is X in units of the absorption length a^^^, and the vertical axis is the pulse 
Area in units of tt. The shaded zone indicates the location of the medium. The solid 
curve is 9a{x)/iT and the dashed curve is 9b{x)/Tr. The top plot shows the 2n Area 
input pulse on the 1 to 3 transition, converted to a 2tt Area output pulse on the 2 
to 3 transition. The bottom curve show the input 1.37r Area pulse first converted 
to a 2tt Area pulse consistent with the Area theorem. The transfer process then 
occurs as predicted by the analytic solution. 



6 Conclusions 

We have presented an analytic solution to the coupled Maxwell-Bloch equa- 
tions that describe a Raman transfer process between two SIT solitons. 
This process can be broken up into three distinct regimes of which only 
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one exhibits the characteristics of a three-level medium. In both asymp- 
totic regimes the process can be considered as a two- level SIT process. We 
also show how to estimate the distance that the pulses will propagate in 
each regime. By knowing the input Areas of the Stokes pulses, one can 
use the Area formulas in Eqns. ^ to determine the propagation distance 
the pump pulse will propagate in regime I, as an SIT soliton. The analytic 
pulse solutions in Eqns. ([5]) then describe the Raman transfer process in 
regime II. Finally, once the transfer process has occurred, we show that the 
Stokes pulse travels as an SIT soliton in regime III for the remainder of its 
propagation through the medium. 
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